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Abstract. The Gradient Echo Memory (GEM) scheme has potential to be a suitable 
5v protocol for storage and retrieval of optical quantum information. In this paper, we 

review the properties of the A- GEM method that stores information in the ground 
states of three-level atomic ensembles via Raman coupling. The scheme is versatile in 
that it can store and re-sequence multiple pulses of light. To date, this scheme has been 
implemented using warm rubidium gas cells. There are different phenomena that can 
^ influence the performance of these atomic systems. We investigate the impact of atomic 

2 motion and four-wave mixing and present experiments that show how parasitic four- 

wave mixing can be mitigated. We also use the memory to demonstrate preservation 
of pulse shape and the backward retrieval of pulses. 
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1. Introduction 

Coherent storage and manipulation of Hght has attracted a lot of attention in the last 
decade due to its application in optical quantum communication p.J and quantum 
computation [2j. There exist various light storage protocols based on coherent 
interaction of light with atomic ensembles. These include elect romagnetically induced 
transparency (EIT) [3], the atomic frequency comb (AFC) |4l[5l[^, four- wave mixing 
(FWM) [46], Raman absorption [7], Faraday interaction [8j, controlled reversible 
inhomogenous broadening (CRIB) [21 HO], and modified photon echo [llj. 

The principle behind CRIB is to construct a reversible absorption process. After 
absorption in an inhomogenously broadened atomic medium, some mechanism can be 
used to invert the detunings of the individual absorbers that are spatially distributed in 
the medium. This reversal in detuning gives rise to a photon echo and thus the retrieval 
of stored optical information. CRIB was proposed in gas cells [13j and solids [Hllll], and 
first demonstrated in 2006 using a cryostatic ensemble of two-level rare-earth ions [15j. 
It was then realized in 2008 that the application of a detuning gradient longitudinally 
along the length of the storage medium allowed for high recall efficiencies in the forward 
direction [ll6j by preventing reabsorption of the light. This overcame a limit of 54% 
recall efl&ciency that had previously been determined \T7]. The use of a longitudinal 
gradient to control the broadening of the atomic ensemble is referred to as longitudinal- 
CRIB or gradient echo memory (GEM). In 2010, ensembles of cryostatic rare-earth ions 
were used in a GEM system to demonstrate the first unconditional quantum memory 
with 69% recall efficiency [18j. 

The focus of this paper is A-GEM, which works using a three-level 'A' atomic 
transition. It was initially implemented in a warm Rb vapor cell in 2008 [19]. 
This scheme has been used as a random access memory for light pulses, [20], has a 
demonstrated efficiency as high as 87% [21j, and provides noiseless storage of quantum 
states [22J. 

In this paper we review the basics of the A-GEM scheme and its description as a 
polariton in the spatial Fourier domain. We discuss the particular issues that relate 
to the gaseous rubidium vapor that is the basis of our experiments, in particular the 
effects of diffusion. We present experiments showing the forward and backward retrieval 
of stored pulses as well as pulse shape preservation. We then model the impact of four- 
wave mixing in our system, which is a potential source of noise for stored quantum 
states. Our work shows, both through models and experiments, that we can mitigate 
the effects of four-wave mixing on the A-GEM scheme. 

2. Light storage using GEM 

Consider an ensemble of two-level atoms, as shown in Fig.[T] Due to the atomic frequency 
gradient, each frequency of the input optical pulse (probe) will absorbed by resonant 
atoms at different points along the length of the ensemble. The atomic polarization 
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along the z direction is thus proportional to the Fourier spectrum of the input probe 
light. To release the stored light, the gradient t] is simply inverted at time r and the 
optical field is regenerated as a photon echo at time 2r. The echo pulse emerges in the 
forward direction but shape-reversed. 

The GEM scheme relies on inhomogeneous broadening being introduced as a linear 
atomic frequency gradient, (^(^,t) = '^(t)^, along the length of the storage medium, 
< z < L, where r\{i) is the slope of the gradient that can change in time. In theory, 
GEM can reach 100% efficiency in the limit of high optical depth (OD) p!6l [23] . 

Two-level GEM has been demonstrated using a cryogenically cooled solid state 
system. A storage efficiency of 69% for coherent states has been demonstrated [18j. 
Furthermore, it was shown that the noise introduced by the memory is negligible for very 
weak coherent states. Using optical transitions, however, means that the storage time is 
limited by the excited state linewidth. Below, we introduce the A-GEM scheme that, in 
principle, offers high efficiency, long storage time, and extra control over manipulation 
of the optical information. 
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Figure 1. Schematic representation of a two- level GEM. An intensity modulated light 
pulse, (a), is incident on an atomic medium, (b), with an atomic frequency spectrum 
{r\z). The bandwidth of the atomic frequency gradient covers the input pulse spectrum 
so that each frequency component is stored at different position along the memory. To 
release the pulse, 77 is switched to —77, (c), after a time r. The stored light emerges as 
an echo, (d), in the forward direction at t = 2r. 



Let us now consider the three-level structure depicted in Fig. [2] (a) where a strong 
coupling field and a weak probe field, with Rabi frequencies of Vt and gE^ respectively, 
interact with a three level atom with one photon detuning A and two-photon detuning b. 
In the weak probe regime, we can assume that once all the population is pumped to the 
state |1) (cTii "r^ 1) it stays there during the interaction time. The Heisenberg/Maxwell 
equations can be written as 
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(Ti3 = - (7 + 7o/2 + iA)ai3 + ig£ + iVtai2 + F13 

^12 = - (70 + ^5(^, t))cfi2 + if^*0-i3 + F12 (1) 

where A/" = gN /c is the effective hnear atomic density. The Langevin operators 
Fi3 and F12 account for noise arising from spontaneous emission 7 and ground state 
decoherence 70 respectively. It has been demonstrated that in the hnear weak-probe 
regime, this system does not show any excess noise beyond that required the preserve 
commutation relations [2^. In other words, the memory can only introduce vacuum 
noise associated with loss - just like a beamsplitter. This means we need not treat the 
noise terms with a full quantum model. Rather, we can run a semi-classical model of 
the system using the equations without the Langevin noise, and then add the vacuum 
associated with the modeled efficiency. In future equations we will omit the Langevin 
terms and the hats from the operators and consider semi-classical c-number equations 
that we can solve numerically. 

Using the steady state solution of Eqs. [T| the susceptibility of the medium can be 
written [25j as 



^ . ^,J2^ + 270(0^ + 707) , (4,5(0^ - 4M) - 4A7g 
^ '|Q2 + (^ + 2^A)(7o + 2^5)|2^ |1^2 + (^ + 2^A)(7o + 2^5)|2. ^ ' 

The normalized transmission of the probe field as well as the imaginary part of the 
susceptibility is depicted in Fig. [2] as a function of the two-photon detuning, 5. At 
large ODs, the absorption of the unbroadened line (trace (i)) is significant and applying 
gradient-broadening of 200 kHz to the narrow absorption feature does not drastically 
alter the maximum absorption (trace (ii)). The real part of the susceptibility, shown in 
Fig. [2](c), shows that there is substantial dispersion of probe field across the absorption 
feature. 

We now simplify the equations further by performing an adiabatic elimination of 
the excited state [26j and assuming ddi'^/dt <C 7, or equivalently 1/T <C 7, where T is 
the fastest time-scale of the evolutions. We also assume a large detuning compared to 
the spontaneous emission rate (A ;» 7). Assuming the coupling beam to be real and 
combining Eqs. [T] we reach 

c^i2 = (-70 + ^5(^, t) - ^^)c^i2 - i^£ (3) 

The term f]^/A is the ac-Stark frequency shift caused by the coupling field, which can 
be cancelled by adjusting the coupling field frequency. Performing the transformation 
£ -^ ^g2fi'A/c/At ^^^ ^ _^ Q^-igNc/At ^^ ^^^ rcmovc the first term on the right hand side 

of the Maxwell equation. We also perform the transformation z^ = z + ct^ i.e. moving 
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Figure 2. (a) Atomic level structure of a A-type transition with two laser fields 
(coupling and probe fields) with Rabi frequency of Q and gS^ respectively, where 
Q ^ g£. The atomic spin decoherence rate is 70, one-photon and two-photon detunings 
are A and 5^ respectively. The excited state linewidth is 7. In the limit of large 
detuning the excited state can be adiabatically eliminated and the A-atom (atom on 
the left) is equivalent to a quasi two- level atom, as shown on the right, (b) Normalized 
transmission of the probe signal as a function of two-photon detuning for unbroadened 
(i) and broadened (ii) Raman lines, (c) The real part of susceptibility corresponding 
to the broadened Raman line. Parameters used to calculate susceptibility are: 1^ = 20 
MHz, 70 = 5 kHz, A = 2 GHz, and r]L = 200 kHz for the broadened hue. 
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to a frame with speed of Hght, to reach 

o"i2 = - (To + i^iz' - ct, t))ai2 - i^S (5) 

The equations for two-level GEM are [27] 

cri2 = - (712 + iS{z, t))ai2 - igS (7) 

^^ = iXai2 (8) 

which are formally equivalent to our Raman scheme if we set Af -^ A/'f^/A, g -^ g^/^ 
and 7 ^ 7o. Therefore, one can see that the GEM storage mechanism described above 
for two-level atoms is also applicable to the three-level systems in certain regimes. The 
two levels are now the long-lived ground states of the three-level system. 

3. Polaritonic description of GEM 

It is possible to associate a quasi-particle (polariton) with the light field and spin wave 
inside the GEM system, similar to the EIT dark-state polariton [3j. The evolution of 
the polaritonic mode inside the GEM system provides a clear intuitive picture of the 
system. Neglecting the decoherence and taking the spatial Fourier transform of Eqs. [7| 
and HI we arrive at 

d d 

--cri2(t, k) = rj—auit, k) + ig£{t, k) (9) 

k8{t,k) =Afai2{t,k) (10) 

and therefore 

Here k is the spatial frequency component of the joint system. A similar equation can 
be written for the electric field, £. Thus, a single mode polariton-like operator in time 
and k space for a two-level atom is obtained [28], ip{k^t) = k£{k^t) + A/'ai2(fc, t), which 
has the following equation of motion. 

Another solution for the above equation of motion as '0(/c, t) = k£{k^ t) —Afai2{k^ t). 
Using EqfTol it is evident that the anti-symmetric mode will never be excited. The 
normal mode equation of motion (Eq. [l2]) indicates that ip{k^t) propagates along the 
fc-axis with a speed defined by the slope of the gradient. As the polariton reaches higher 
k values the electric field amplitude decreases. As each frequency component of light 
propagates through the medium and gets absorbed by resonant atoms, the dispersion 
due to the neighboring atoms, having only slightly different resonant frequencies, will 
affect the propagation and absorption of the light. This is, in fact, the origin of the 
photonic component of the polariton and also the phase shift of the echo investigated 
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in Ref. [12j. The intensity of the atomic excitation, created after the input hght enters 
the medium, remains unchanged during the storage time. Fhpping the gradient (i.e. 
T] ^ —rf) time reverses the absorption process so that when the polariton reaches fc = 
(the phase matching condition) a photon echo emerges from the ensemble in the forward 
direction. We note here that at any time, the system is described by a normal mode 
with different amplitudes for optical and atomic fields, therefore it makes more sense to 
picture the GEM storage as continuum of evolving polaritons rather than just a single 
polariton. 

In a three-level atomic ensemble, in the far-detuned and adiabatic regime, using 
Eqs. ([5]) and ([g]) a similar polaritonic mode can be described as ip{k^t) = kS{k^t) + 
A/'Vi2(A;, t), where A/"' = J^Q/A. The second term in this equation necessarily goes 
to zero if the coupling field is switched oS during storage. Therefore, as it is still the 
case that k£ = A/'Vi2, the slowly propagating electric field vanishes and the polariton 
becomes purely atomic. 

The real parts of the electric field and atomic polarization, for a modulated input 
pulse, are plotted in Fig. [3] (a) and (b), respectively, obtained using an XMDS [29j 
numerical simulation of Eqs. ^ and ([s]). It can be seen that the light is nearly 
stopped at the centre of the ensemble after the pulse enters the medium and its intensity 
gradually decreases. The atomic spin wave is quickly generated after the pulse enters 
the memory. We also see that the spatial structure of the atomic and electric fields 
becomes increasingly fine as time progresses before gradient switching. 

The spatial cross section of the atomic polarization at any time during the storage 
is the Fourier spectrum of the input pulse. This explicitly demonstrates the frequency- 
encoding nature of GEM. Conversely, any cross-section of the polariton ip{k^t) along 
the k axis shows the temporal profile of the pulse [28] as a second Fourier transform 
into fc-space returns the original pulse shape. This is shown in Fig. [3] (c) where the 
polaritonic normal mode of the modulated pulse is plotted. 

It is quite interesting to compare the polariton pictures in EIT and GEM 
systems [28j . The EIT polariton propagates through the medium {in z — t plane) while 
the control field is on and becomes stationary when control field is switched off. Thus, the 
temporal modes compressed within the atomic sample and, by integrating the polariton 
amplitude along the propagation axis, one can construct the temporal shape of the light 
field. The GEM polariton, on the other hand evolves in the spatial frequency domain 
and its evolution can be controlled via the external field, as described above. 

The group velocity of the light field propagating inside the GEM medium is given 
by Vg = gAf/k'^. The further away the polariton is from A; = 0, the less intense the 
electric field and the smaller the group velocity becomes. In the A-GEM scheme, the 
group velocity and light amplitude can be controlled by both the coupling field and the 
gradient field. This extra control over the information provides considerable fiexibility 
for data manipulation. 
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3.1. Steering of the GEM polariton 

The dynamic characteristics of the GEM polariton aUow one to precisely control and 
manipulate the state of the system in time. This control is even more versatile in A- 
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Figure 3. Numerical simulation of the storage of an intensity modulated pulse. 
(a) The real part of the electric field m. z — t plane. Inset shows the temporal shape 
of the pulse at the input and at the output. The electric field decreases in strength 
until the atomic frequency gradient is switched dX t — 7b \is. It then increases and 
eventually an echo is emitted, (b) The real part of the atomic polarization in the 
z — t plane. A cross-section of the atomic polarization along the z-axis reveals the 
Fourier spectrum of the input probe light. The inset shows the absolute value of the 
polarization at the indicated position, (c) Intensity of the polaritonic excitations in 
the k — t plane, /c = indicates that the rephasing of the atomic spin is completed 
to produce a coherent emission. A cross section of the polariton along k reveals the 
temporal shape (amphtude) of the stored pulse (as shown in the inset). Parameters 
used for these simulations are: gM /r] = 1.74, gL = I67, and 70 = 0. 
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GEM due to the contribution of the couphng field and detuning gradient in determining 
the state of the system. 

3.1.1. Atomic detuning. The detuning gradient causes the different frequency 
components of the hght to be absorbed spatiaUy along the atomic memory. In the 
polaritonic picture, the phase matching condition is equivalent to A; = 0. After the 
light is stored, r]{z) can be switched to any arbitrary shape in time. The photon echo, 
however, is only emitted when the atoms are rephased and coupling field is on. The 
shape of the frequency detuning is chosen to be monotonic in order to avoid any re- 
absorption. One can use other forms of frequency detuning to avoid reabsorption, but 
a linear detuning is practically easier to implement and control. 

The shape of the detuning can be engineered in such a way as to optimize storage 
of a light pulse with a particular frequency spectrum. The size of the frequency gradient 
determines the bandwidth of the memory and, by applying a larger gradient field, 
the storage bandwidth increases. The number of well separated pulses, each with a 
bandwidth smaller than the memory bandwidth, that can be simultaneously stored is 
limited by the memory lifetime. For a simple linear gradient, the efficiency of GEM 
is inversely proportional to the size of the gradient. There will, therefore, be a trade- 
off between the bandwidth and storage efficiency. For efficient storage of broadband 
information, large ODs are imperative. 

The frequency gradient can be switched to higher or lower values or even to zero 
during the storage time without affecting the storage mechanism. In Fig. |4] we show 
numerical simulations for the evolution of the normal mode in the k — t plane where tj 
is switched to different values. In Fig. [i] (a), a recall gradient four times steeper than 
the input gradient introduced and therefore the atomic field evolves faster towards the 
origin. Shortly after, before the polariton reaches /c = 0, ry is switched to zero to stop 
the evolution and freeze the polariton at a constant k value. The system can then be 
described as a single polariton. Eventually, with tj = — ryo, light is emitted when atomic 
excitations are phase-matched at k = 0. 

3.1.2. Tunable atom-light coupling. The role of the classical coupling field in A-GEM 
is to couple the weak probe light via a Raman transition to the ground states. The 
coupling field is also required at the reading stage when the photon echo is emitted to 
transfer the information from the atomic spin back to the light field. In an ideal system, 
the dynamics of the coupling field during the storage time would not affect the storage 
process, but practically it is beneficial to turn off the coupling field during the storage 
time to reduce scattering and therefore loss. The scattering process is further discussed 
in the experiment section. 

The coupling field intensity determines the effective atom-light coupling and 
therefore the OD. Decreasing the coupling field intensity reduces the storage efl&ciency at 
the writing stage. At the reading stage this can result in partial retrieval; the intensity 
of the coupling field when the normal mode reaches fc = determines how much stored 
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Figure 4. Numerical simulation of GEM showing absolute value of the polaritonic 
excitation in k — t plane, (a) Multiple switches of the atomic gradient field. The 
switching protocol of the gradient detuning is depicted on top where r] is switched to 
T] = —47^0, then r? = 0, r] = Ar] and finally to r] = —r]^. (b) Arbitrary probe retrieval. 
The coupling field is turned on during the writing stage and only for a short period 
during the reading stage when the condition /c = for the middle pulse is satisfied. 
This guarantees that only the second pulse is recalled. The top part of the figure shows 
the switching protocols for coupling field intensity , /c, and detuning gradient, r]. The 
inset shows the temporal profiles of the input and echo pulses. 



excitation is converted to the electric field and how much is left inside the memory. 

Fig. |4](b) shows a simulation of the normal mode evolution in the k — t plane in 
the A-GEM system. The coupling field intensity is switched to zero shortly after the 
three light pulses enter the medium and their information is transferred into the atomic 
excitations. After the gradient is flipped and the excitations reach fc = 0, the coupling 
field is turned back on only for the time window t = 16 to 187"^ in order to couple out 
the second pulse. In this case, the information imprinted into the atomic ensemble from 
the two other pulses will remain inside the memory. 



4. Time sequencing and manipulation of pulses 

Using the polariton description introduced above, we show here how optical information 
can be coherently manipulated in an arbitrary manner. We discuss the possibility 
of coherent pulse sequencing, with first-in-last-out (FILO) , first-in- first-out (FIFO), 
arbitrary recall and backward retrieval of information using A-GEM [20j. 
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4.1. First-in-last- out storage 

One of the intrinsic properties of GEM is that the photon echo is a time-reversed copy 
of the input pulse. When the optical field frequencies are stored spatially along the 
z-axis, fiipping the gradient triggers the time-reversal process so that the pulse sequence 
is reversed. 

Fig. [5] (a) shows numerical simulation of the polariton for first-in-last-out (FILO) 
storage of four input pulses, depicted in k — t plane. The first pulse to enter the system 
returns to A; = last, and is thus re-emitted last. 

4.2. First-in- first- out storage 

Shape preserving storage in GEM is possible if the absorption is time-reversed twice. 
The first pulse in would then be the first pulse out. Using a three-level system, one 
is free to turn off the coupling beam. In this case, although the atomic polarisation 
is completely rephased when the normal mode reaches A; = 0, no light can be emitted 
because the atom-light coupling is turned to zero. This is seen in Fig. [5](b), which 
shows |'0(t, /c)p for the switching scheme shown on the top part of the figure. With the 
coupling beam off, the normal mode passes straight through k = to negative k values. 
We can then switch the frequency gradient once again to time-reverse the atomic spins 
one more time. Now with the coupling field back on, the normal mode is converted into 
a photon echo at A; = without pulse shape reversal. In this way, we can construct a 
first-in-first-out (FIFO) memory. 

FIFO storage can also be achieved in a single memory operation based on two-level 
GEM, albeit with some loss. This can be done by switching the gradient to a very steep 
value and opposite sign after storage as shown in Fig. [5] (c). This will effectively reduce 
the OD of the system at the first reading stage and therefore most of the excitation 
will remain inside the memory j30j. The left-over excitations, that contain most of the 
information, can be recalled later by applying a gradient similar to the one at the writing 
stage. 

4.2.1. Multiple retrievals In the two-level GEM scheme, for low OD regimes, the 
writing stage becomes inefficient as light will leak through the atomic ensemble without 
being absorbed. The fraction of the light that is absorbed is 1 — exp(— 2/37r), where 
(5 = gAf/rj. Similarly, at the read stage, low OD means that not all of the atomic 
polarization will be converted back to optical field. A fraction exp(— 2;57r) of the atomic 
excitation remains trapped in the memory [30j. Provided that the coherence time of 
the atoms is long enough, part of the left-over excitation can be converted back to a 
light field by fiipping the field gradient multiple times. In fact, multiple switching of the 
gradient allows atomic excitations to travel back and forth in k space, and each time they 
pass through fc = 0, a fraction of excitation can contribute to photon echo emission until, 
eventually, all the left-over excitation becomes depleted. Disregarding decoherence, the 
fraction of the input light released in the n^^ echo is [1 — exp(— 2/37r)]^ exp[— 2/57r(n — 1)]. 
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Figure 5. Numerical simulation of absolute value of the polaritonic excitations 
showing, (a) first-in-last-out and (b) first-in-first-out in a three-level system of four 
input pulses m. k — t plane, (c) FIFO memory in a two-level system. The original 
gradient is switched to a much steeper negative gradient to reduce the effective OD of 
the memory. The two pulses are retrieved in FIFO manner after switching the gradient 
once more to its original slope. In this case, around 7% of the light is retrieved as FILO 
at the first readout stage, (d) Storage of a train of four pulses in 1,2,3,4 order and 
recalled in 3,2,1,4 order. The switching scheme of the coupling field intensity and 
gradient detuning is plotted on top of each part. Insets show the temporal shape of 
input (blue) and echo (red) pulses. 



In the case of the A-GEM scheme, however, the presence of the coupling field 
provides another degree of freedom to control the amount of atomic excitation converted 
to the light after each gradient flip. By leaving the coupling field off when fc = 0, the 
emission is suppressed and the atoms hold on to the information. Partial retrieval can 
be achieved by lowering the coupling field power at the readout. This ability can be 
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used to create a time-delayed beam splitter or polariton-polariton interferometer J3T] . 

4.3. Arbitrary retrieval 

By controlling the frequency gradient and coupling field intensity in time, it is possible 
to construct a system that can recall the pulses in any arbitrary order. In the case of 
multi-pulse storage, one can decide to turn on the coupling field only during the chosen 
time windows, when atomic excitations are at fc = 0, to recall only desired bins of the 
stored information and leave the rest inside the memory. A decoherence-free model of 
the on-demand retrieval is shown in Fig. [5](d), where 4 pulses are stored in the memory, 
and after the first field switch, pulses 3 and 2 are recalled by turning the coupling field 
on, only during the time window that these two pulses cross k = 0. The other two pulses 
are recalled later after the second field switch. One can choose any other combination 
in the reading stage by controlling the detuning gradient and also the coupling field 
switching. 

4^4' Backward propagating retrieval of echo 

Photon echoes created using two-level GEM always co-propagate with the input pulse. 
The A-GEM is capable of recalling stored information in both backward and forward 
directions [32j. By turning off the coupling field co-propagating with the signal field, 
the slow light vanishes, and a collective spin excitation is created. At the reading stage, 
when the coupling field is turned back on, the signal pulse is recreated, propagating 
along the quantization axis in the same direction as of the coupling field. If a coupling 
field counter-propagating with the initial signal field is applied at the read out stage a 
photon echo will be emitted in the backward direction. We demonstrate this below in 
the experimental section. It is also possible to create stationary light inside the memory 
by means of counter-propagating coupling fields, similar to EIT [33l [34]. This matter 
will be left for the future publications. 

5. Implementation in warm vapor 

We have experimentally implemented A-GEM using a ^^Rb vapor cell at 70-80°C At 
such high temperatures the effects of Doppler broadening, diffusion, and collisions can 
impact the storage mechanism. Before presenting experimental data, we will discuss 
some aspects of warm vapor that impact on the implementation of GEM in Rb vapor. 

5.1. Angular dependency of absorption and transmission lines 

The mode-matching and crossing angle between the coupling and the probe beams are 
crucial for performance of both EIT- and Raman-based memories. To study this effect 
we consider the angular configuration of optical and atomic beams shown in Fig. [6] where 
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Figure 6. Schematic diagram showing crossing angle between optical and atomic 
beams. The reference coordinates are v\\ and v_\_ representing quantization axis and 
axis perpendicular to it, respectively, m and Va are the mass and velocity of the atom, 
kp is the wave- vector of the probe light and kc is the wave- vector of the coupling light. 



the angle between the control and probe is shown as 9 and the atomic velocity makes 
an angle with the quantization axis. 

In a vapor cell, if the crossing angle between the beams is non-zero, the one-photon 
and two-photon detunings need to be modified to include extra frequency shifts [35] in 
order to describe the physics of the system more accurately. Change in the one-photon 
and two-photon detuning of the Raman line are given by 



A^ A 

5 ^ 5 + (kp 



kp.v 



kc).v 



(13) 



where 



kp.v 



rhpVp 



= - kpVa cos{9/2 + (/)) 
(kp - ke).v = 2kpV^ sm{9/2) (14) 

and bold letters represent vectors. The velocity, Vp^ is the component of atomic velocity 
along the probe beam and v± is the component of atomic velocity perpendicular to the 
quantization axis. Taking these frequency shifts into account, the imaginary part of the 
susceptibility (see Eq. [2]) can be rewritten as 

(8(5- Vp)'7 + 27o(^' + 7o7) 



Im[x\ 



(15) 



\n^ + (7 + 2z(A + kpv^ sin(e/2)))(7o + 2i{5 - Vp))P 

At temperatures higher than room temperature, the atomic velocity in each 

direction is defined by the Maxwell-Boltzman distribution. We assume that the velocity 

is isotropic in the vapor cell, so the distribution of v^ is the same as the distribution of 
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v± and Vp. The absorption coefficient of the medium is then given by 

/oo roo 
/ Im[x\{vp,v±)N{vp)N{v±)dVpdv±, (16) 

-co J —oo 

where N{v) is the Maxwell-Boltzman distribution function. This integral can be 
numericaUy solved at T > to find the absorption coefficient and the FuU-Width- 
at-Half-Maximum (FWHM) of Lorentzian-shaped Raman lines. The Raman linewith 
shows a linear relationship with the crossing angle as shown in Fig. [7] (a). The absorption 
coefficient for two different values of f^, as a function of the crossing angle, is plotted in 
Fig. [t] (b). The crossing angle corresponding to the maximum absorption of the Raman 
line approaches zero at A/Q ^ 1. 

5.2. Doppler broadening 

In order to describe the properties of warm vapor cells used in our experiment, we need 
to consider Doppler effects caused by the thermal motion of the atoms. The velocity 
spread of the contributing velocity classes can be roughly determined by taking into 
account the optical linewidth of the excited stated denoted by 7. The velocity spread 
leads to a Doppler broadening of the observed Raman line in the laboratory frame. This 
width at T c^ 70°C is around F^ = 500 MHz. 

We follow a similar method as the one described in Ref. [36j and solve the integrals 
numerically to calculate the FWHM of the Raman line for different coupling field 
Rabi frequencies. The results of numerical calculation of the Raman linewidth for two 
different detunings as a function of Q are depicted in Fig. [t] (c). We note here that the 
far-detuned two-photon absorption, i.e. A ^ T^^ would be insensitive to the Doppler 
broadening of the excited state as long as the two fields are co-propagating. 

5.3. Buffer gas in the vapor cell 

In most of the light storage experiments that have been performed on warm Rb vapor 
cell, inert buffer gases such as Ne, Kr and He have been used to reduce the mean free path 
of the Rb atoms and therefore increase the time of ffight of atoms inside the laser beams. 
The collision of Rb with buffer gas is elastic and therefore the phase information is better 
preserved, i.e. the ground state atomic spin remains largely unchanged. Collisions with 
buffer gas may, however, shuffle population in the excited state and can cause collisional 
broadening [37] that is given by 

where p is the gas pressure, T is the temperature and v is the average velocity of atoms 
and cT/e is the kinetic cross section. At room temperatures, the collisional broadening for 
Rb-Kr is Ff^^-^^ = 17.1 MHz/Torr and for Rb-Ne is Ff^^"^" = 9.84 MHz/Torr [38j. We 
use 0.5-1 Torr of Krypton as buffer gas in our experiment due to its similar mass to Rb. 
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Figure 7. (a) FWHM of Lorentzian Raman absorption line calculated using Eq. 16 
where Q = I.O7. The solid line is to guide the eye and fluctuation of the points is due 
to lack of precision in the integration, (b) Normalized absorption coefficient calculated 



for detuned light by numerical integration of Eq. 16 for different crossing angle of 
coupling and probe beams. Red (i) and blue (ii) data show the absorption coefficient 
for ft = I.O7 and Q = 2.O7 respectively. The dashed lines are guides only, (c) FWHM 
of the Raman line plotted using numerical integration for ^ = 0, (i) A = IOO7 and (ii) 
A = 3OO7. Doppler effect has been included in the model and other parameters used 
to plot (a), (b) and (c) are: 70 = 0, T = 350K. 
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5.4' Atomic diffusion in a gas cell 

At low buffer gas pressures, a major relaxation process appears due to spatial diffusion 
out of the laser beams. In our case, we use a cylindrical geometry with transverse 
diameter of a few millimeters. We are interested in how long it takes for an atom to 
diffuse out of the transverse cross-section (i.e. the detection mode volume). 

In two dimensions, the standard deviation in position a^ (defined such that 
approximately 32% of atoms have moved a distance greater than a^) is given by [39j 



ar = 2V2Di (18) 

where L^ is a diffusion coefficient with units of cm^/s that is defined as 

D = Do{^) (19) 

Here Pq = 760 Torr, P is the buffer-gas pressure in Torr, and t is the time. The diffusion 
constant Dq for Rb-Kr is calculated at temperature of 45° to be about 0.16 cm^/sec [40j. 
The diffusion constant above room temperature is scaled with T^/^. 

Longitudinal diffusion could also play a role, especially considering the longitudinal 
frequency gradient. If an atom moves longitudinally by a distance 5z during the storage 
time, it will experience a frequency shift of rjSz at read out. This means that longitudinal 
diffusion will lead to random frequency shifts of the atoms if the applied gradient is 
sufficiently large. The effect of longitudinal motion in our experiment is negligible 
because r]6z is small compared to the unbroadened Raman line. 

6. Experimental setup and results 

The results presented in this section were obtained using a ^^Rb mixed with 1 Torr 
of Kr buffer gas in a warm cell with length of 7.5 cm and diameter of 2.5 cm. The 
experiment setup is shown in Fig. [8} The probe and coupling fields were collimated 
with 3 mm and 5 mm diameter, respectively, and had orthogonal linear polarization. A 
second counter-propagating coupling field was used to retrieve the echo in the backward 
direction. The Raman detuning was around 1 GHz and the power of the probe and the 
coupling fields were typically around 10 |llW and 200 mW, respectively. To generate 
the linear field gradient we used two specially designed solenoids and, by switching the 
current between them, the sign of the gradient is switched [2T]. We used two layers of 
/i-metal shielding to reduce the background magnetic field. 

The coupling field can drive the population out of the level |1) or |2) and cause 
spontaneous Raman scattering (SRS). This can reduce the coherence time of the 
memory. The loss due to the scattering is given by Tgcatt = 7a/1 + f^^/(7^ + A^). 
At large detunings, the scattering rate follows the laser power. At low coupling field 
powers the FWHM of Raman line is the ground state linewidth 70. We observed that 
the switching off the coupling field during the storage time increases the coherence time 
by almost an order of magnitude ^21j . 
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Figure 8. Schematic experimental setup for forward and backward retrieval of 
echo signal from the memory. The probe field is shifted by 6.8 GHz using a fiber 
coupled electro-optic modulator (FC-EOM) to match the frequency splitting between 
the hyperfine ground states of ^'^Rb. A cavity was then used to separate the 
carrier and -6.8 GHz sideband from the probe field after the FC-EOM. BS: beam 
splitter, PBS: polarized beamsplitter, AOM: acousto-optic modulator, 8^^: backward- 
propagating probe field, f/^: forward-propagating probe field, 8^^: backward- 
propagating coupling field, f/^: forward-propagating coupling field, ^6//: detectors 
for the forward and backward-propagating modes. 



6.1. Results 

6.1.1. Shape mirroring One of the intrinsic properties of GEM is that the recalled 
photon echo has the same shape as the input signal but time reversed. FigJ9] (a) shows 
the experimental proof of the shape mirroring for a double Gaussian and ramp shape 
input pulse, respectively. Output echoes are indeed a mirror image copy of the input 
pulses due to the time-reversed nature of the process. 



6.1.2. Backward retrieval of echo To experimentally demonstrate backward retrieval 
of the echo, we used a coupling field propagating in opposite direction with respect 
to the input pulse. Initially, the input pulse is stored with a co-propagating coupling 
field. At the read out stage the photon echo was retrieved either in forward or backward 
direction by means of forward or backward coupling field after the gradient flip. The 
results are shown in Fig. |9](b) where photon echoes retrieved in the forward (trace (ii)) 
or backward (trace (iii)) directions. 



7. FWM and Raman absorption in a dense atomic sample 

The FWM phenomena in EIT systems was studied in Ref. [41] and experimentally 
observed later in a dense atomic sample [121 US] • It was shown that FWM in EIT can 
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Figure 9. (a) The shape mirrored photon-echo for a train of two pulses with different 
amplitude ratio. Inset shows a ramp-shape input (i) and shape mirrored photon echo 
(ii). Top section shows the switching protocl for the coupling fields, (b) Forward 
and backward retrieval of echo, (i): The input pulse, (ii): Forward propagating 
signal light at the output of the memory when the coupling and signal beams are 
co-propagating, (iii): Backward propagating light at the output of the memory when 
a coupling field, counter-propagating with respect to input hght, is applied at the 
retrieval stage. Backward propagating light measured using a 50/50 beam splitter 
before the memory, for this reason trace (iii) is amplified by a factor of 2. Top section 
shows the switching protocol for the forward (blue) and backward (red) propagating 
coupling fields. Due to the low OD, part of the input light leaks through the memory 
without absorption that can be seen in both (a) and (b). 



limit the storage efficiency at high optical depth [42j. In this case, the conversion of 
an original pulse from the signal to Stokes channel may reduce the read-out efficiency 
[44]. Under certain conditions, however, FWM may lead to gain for both the signal 
and Stokes fields, which could compensate for any optical losses [45j, albeit with added 
noise. 

With high OD, FWM phenomena can also be observed in a GEM. We consider the 
level scheme shown in Fig. [lO] to derive the equations of motion. Assuming a linearly 
varying two-photon detuning S{z) = rjz and a far-detuned Raman transition, A :» 7 we 
can arrive at the following equation for probe {£) and Stokes {£s) fields at the centre 
frequency 

dz \ S:{z) J V «21 a22 J V ^siz) J 



(20) 



N 



oo 



7c(i]2 + rro) 
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an = igT(j 



ai2 — a2i 



a22 = - ir 



g'^a 



/0/2 



A'2 



(21) 



and 



To = (70 + i5) 
r =(7 + «(A 



5)) 



(22) 



We solve these equations numerically in different regimes to understand the contribution 
of the FWM process in a Raman absorptive medium. 



Figure [TT] shows the calculated transmission of the probe and Stokes fields for two 
different values of resonant 0D= g^NLj^c and two different values of the initial Stokes 
field amplitude, parameterized by r = Es[z = 0)/S{z = 0). In our experiment, a non- 
zero Stokes seed is a possibility. The probe field is formed by one sideband of a phase 
modulation that is filtered by one or more resonant cavities. If the leakage of these 
cavities is too high, the other sideband of the phase modulation can form a Stokes seed. 
It can be seen that for large OD, probe amplification occurs around the edges of the 
Raman line. The results show that 4WM amplification of the probe can be suppressed 
by minimizing the seed light of the Stokes field (making r small) and by using an atomic 
system where Am ^ 7. The amplitude of the Stokes field, as well as the gain, is also 
substantially reduced at the point of maximum probe absorption. The Stokes field is 
also absorbed around two-photon resonance. FWM inside a vapor cell at 140°C has 
been experimentally investigated in a detuned double- A configuration (with A = I5O7) 
[46] . where the maximum gain was observed 20 MHz away from the Raman resonance. 




I1> 



Figure 10. Atomic level diagram for degenerate FWM between coupling and probe 
fields. The Rabi frequency of the coupling field for different transitions are shown by Q 
and Q' due to different dipole moments. The Rabi frequency of the probe and Stocks 
fields are shown by gS and g'Ss^ respectively. 
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Figure 11. Transmitted amplitude of, (a)-(a') probe, (b)-(b') Stokes fields 

normalized to the input probe amplitude for two different ODs. A transmission 
greater than 1 indicates gain in the medium. Parameters used for this numerical 



simulations are: Tp = 47, 70 = O.OO27, ft = 87, OD = g NL/jc 



65,280, 



r = Ss{z = 0)/S{z = 0) = 1.0,0.01 and A = 2OO7. In the limit Ahf > 7 we can 
effectively switch off the FWM process. In these calculations, the Doppler effect has 
not been included in our model 



In a dense atomic medium and at Raman resonance, the competition between large 
amplification and large absorption leads to complex dynamics which can result in a 
breakup of the probe pulse. The interference between the probe and Stokes transition 
can be seen in Fig. 11 (a^). 



Using linear polarization instead of circular polarization can significantly enhance 
the FWM process in a Rb vapor cell. Linear polarization is treated by the atoms as 
a superposition of left and right circular polarization and will therefore interact with 
multiple excited state sublevels. This will effectively increase the atom-light coupling 



strength that enhances the non-linear process of FWM. Fig. 12 (a) shows results of 
numerical calculations for transmitted probe and Stokes signals for two different ODs 
assuming co-propagating beams. 

Experimental results showing the broadened Raman lines for different coupling field 
powers and temperatures are shown in Fig. [12] (b). These experimental results have been 
obtained using a 20 cm gas cell containing ^^Rb and 0.5 Torr of Kr buffer gas, while 
the coupling and probe signal fields have orthogonal linear polarization. The probe 
measurement was performed using heterodyne detection. As shown in Fig. [8] the probe 
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field is one of tlie 6.8GHz pliase modulation sidebands generated using a fiber coupled 
electro-optic modulator (FC-EOM). A cavity of finesse ^ 100 was used to select out 
the probe frequency. Nevertheless, a small amount of seed light (from the other phase 
modulation sideband) will leak through and enhance the amplification process. Using 
circular polarization, we can achieve strong absorption of the probe with no measurable 



gain, as shown in Fig. 12 (c). In this case we have also used two optical cavities for the 
probe beam, each with a finesse of about 100 providing suppression of the seed signal 
of more than 10^ times. 

Applying the detuning gradient increases the absorption bandwidth for the probe 
and therefore widens the region around two-photon resonance in which gain is 
substantially suppressed and the memory can operate noiselessly [47j. At regimes of 
large OD applying a gradient does not noticeably reduce the maximum absorption of 
the probe and it causes suppression of the gain and generated Stokes field. We note 
here that, in the experiment configuration, all of the beams are co-propagating together 
and therefore the phase matching condition is satisfied for FWM process. However, 
the Raman gain can be suppressed if the two A transitions destructively interfere. The 
phase of the seeded Stokes field in the experiment is not controlled and we believe that 
the observed gain is largely due to the vacuum seed. In our model, however, we consider 
phase of the seeded Stokes field such that the gain is maximum. 

We calculate the total energy of the Stokes field around two-photon resonance and 
also 0.047 from resonance (where the amplification is maximum) for different optical 
depths in an inhomogenously broadened medium. As it can be seen in Fig. [T3|(a) and 
(b), the energy of the probe and of the Stokes field around the two-photon detuning, 
^{(jUq = 0) = J_^£{uj — u)o)d(jj^ resonance rapidly decreases by increasing the optical 
depth of the medium. However, 0.047 away from the resonance the probe and Stokes 
energy, ^{ujq = O.O47) = j_^£{uj — (jjQ)d(jj^ increases with OD. 

The significant reduction of Stokes field around the centre of the Raman line 
indicates large suppression of the FWM process and therefore lack of amplification 
in this region. The situation is reversed in a dense EIT medium, where the probe 
transmission is enhanced. 

8. Summary 

In summary, we have discussed properties of A-GEM in a warm vapor cell and presented 
experimental results of shape mirrored and backward retrieved echo. We proposed a 
scheme based on counter-propagating coupling fields to observe stationary light eflFect 
in GEM system. We also discussed the eflFect of FWM in GEM storage process and 
interplay between these two processes in a dense atomic medium. 

This research was conducted by the Australian Research Council Centre of 
Excellence for Quantum Computation and Communication Technology (project number 
CE110001027). 
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Figure 12. (a) Transmitted amplitude of the probe and Stokes fields normalized 
to the input probe amplitude. Parameters used for these numerical simulations are: 
70 = 0.0027, 1] = 3 and 5.27, g'^NL/^ic = 470,1550, r = 0.01, A = 2OO7 and 
r]L = O.O87, T = 85^ C. In these calculations, the Doppler effect has not been included 
in the model, (b) Experimental measurements of the probe transmission for different 
temperatures and coupling field power obtained using linear polarization. Absorption 
is maximum at a non-zero two-photon detuning because of the DC offset magnetic 
field, (c ) The heterodyne signal showing the broadened Raman absorption line for the 
same circular polarization of the pump and probe fields with pump power of 350mW. 
The probe power in all cases was ~ 1 |llW. The envelope shown is a guide to the eye. 



Storage and Manipulation of Light Using a Raman Gradient Echo Process24 




500 1000 1500 2000 2500 

OD 



Figure 13. Total (a) probe and (b) Stokes energy within a frequency window of 
25 = 0.047 around zero two-photon detuning (dashed black line) and uq = O.O47 
away from the two-photon resonance (solid red line) calculated for different OD. This 
quantity corresponds to the total probe and Stokes field energies at the output of the 
medium. Both axes are normalized to the input probe energy. Parameters used for 
this numerical simulations are: 70 = O.OO27, A = 2OO7, Q = 87, r = 0.01, r]L = O.O87. 
In these calculations, the Doppler effect has not been included in the model. 
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